= a step input applied at k r th step (43) 

where 



u n-k t f 0 {n = 1,2,..., hi- I) \ 
1 1 (n > J 



(44) 



To assure independence of the inputs, the onsets of the signals should be apart from 
each other by a sufficient number of steps, i.e., ki must be large enough. 
5 Similarly, in another alternate approach, a pulse input based SCI (PSCI) may be used. 

In this approach, one can also apply multiple pulse inputs in a sequential manner: 

r n = S n ~ ki 
1 % 

= a step input applied at fc r th step (45) 



where 



\ 0 (n = all other points) / [q0) 



10 Again, ki should be large enough to ensure independency of the applied signals. 

Yet another alternate embodiment to the SCI/ERA methods discussed above will now 
be described. In applications of discrete-time computational models, there exist two 
conflicting requirements for the incremental time step A/. For numerical convergence, one 
should adopt a sufficiently small Atj. On the other hand, given the highest frequency of 

15 interest, car, the Nyquist criterion requires < ^. Usually, for practical purposes 
Af 2 » 
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For instance, in a structural model that is coupled with a CFD model for aeroelastic 
applications, the highest mode usually has a much lower natural frequency than the highest 
frequency content in the aerodynamic model. If the signals used in the SCI/ERA methods are 
sharp as in the random, step, or pulse inputs, the SCI will excite all the system dynamics and 
5 hence this characteristic will be carried over to the ERA based reduced-order model. As a 
result, the ERA reduced-order model (ROM) may be prone to have a large dimension to span 
the same high frequency range as the original full-order model (FOM), which suggests that 
there may be room for further order reduction in the ROM. 

To perform a second order reduction, one can apply the Frequency-Domain 

10 Karhunen-Loeve (FDKL) method to the ERA ROM defined by matrices, (28)-(31), wherein 
frequency samples of the system within the given frequency range, (-co 0 co c ) are used to 
extract optimal modes, and a new reduced-order model is constructed via the Galerkin's 
approximation, (see, e.g., Kim, T., Discrete-Time Eigen Analysis and Optimal Model 
Reduction for Flutter & Aeroelastic Damping/Frequency Prediction Based On CFL3D-ERA, 

15 B-ADVTECH-LLL-M02-013, BCAG, February 2003, incorporated herein by reference). In 
this embodiment, the optimal or so called KL modes 0*, are the eigenmodes of the 
covariance matrix K: 



K<k = A,0, (47) 

where 

Kij = X(u i )X(u j y T (48) 
u>i = sampling frequencies 

= uj m ] (49) 
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where coi = -co c and co M = co c . X(nJ are the frequency solutions of the ERA ROM 
subjected to the single-composite-input described by (20) and (21) except that it is prescribed 
in the frequency domain. Once the optimal modes are obtained, assume 

X - *P (50) 

where 

* = l<t>l<t>2 -.. 0«J 

Pi 
P2 

PR, 



P = 



(51) 

Pi 

(52) 



Ri is set to be equal to the rank of the covariance matrix which is usually smaller than 

R. 

After inserting (50) into (1) and (2) with the ERA ROM matrices, premultiplying both 
1 0 sides by * T yields a new reduced-order model as 

p n+1 = Aj p n + B x u n (53) 

y n = Cip n + Du" (54) 

where 

Aj = * T A* (55) 

B, = * T B (56) 

Cj ee C * (57) 

The dimension of the new model is (R/xRi). 

Application of various embodiments of methods and systems to representative, large- 
15 scaled CFD models will now be described. Unlike the general system described by 
Equations (1) and (2), an unsteady fluid dynamic system is driven by displacement and 
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velocity of its moving boundary surface simultaneously as they both contribute to the normal 
downwash on the surface. If one considers a statically nonlinear, dynamically linearized flow 
field, the unsteady fluid motion can be described as 



n+l 



= Ax n + B 0 u n + Bid" (58) 
y n = q(Cx n + D 0 u n + D x u n ) (59) 



where 



x = (Lxl) fluid states (60) 
u = (N { x 1) generalized displacements (61) 
li = (N { x 1) generalized velocities (62) 

y = (Ni x 1) generalized aerodynamic forces (63) 
q = dynamic pressure (64) 

It is noted that the above equations progress in non-dimensional time, r s x > rather 
than in the real time / and ( ) is the first derivative with respect to t. In fact, the dependency 

10 of the normal downwash on air speed disappears when the equations are discretized in t, as 
in Equations (58) and (59). The structural degrees of freedom, w, are the generalized 
coordinates associated with structural modes. These modes are used to describe the motion 
of the lifting surface. Two different types of reduced-order fluid dynamic models can be 
obtained depending on how the inputs are treated. If necessary, the FDKL/SCI can be 

1 5 performed for additional reduction. 

In one particular embodiment, a method in accordance with the invention may be 
applied to an aerodynamic ROM with displacement and velocity inputs. In this embodiment, 
one can treat u n and u n separately as independent inputs. Thus, for the pulse inputs 
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{5 ,„. ( u:%} (65) 

for i = 1,2,...,N U we obtain the corresponding responses y?,yj at the first two 
time steps. Let us define 

Y° = [y?y5...y?v i y^ + iy 0 Ar i+2 ---yk] ^ 
Y 1 = [y!y^...y] Vj y I V i+ iyk + 2---y2N i ] ( 67 > 

where the first N { samples are due to the pulses in u n and the next N { ones are due 
to the pulses in u". Next, we prepare the following inputs, 

»=1 i=l 

d n sci = Edoi-f + Edi,rR l44 (69) 
»=i fei 

10 Subject to the SCI, we sample the system response y" and get 
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2Ni 

y n+1 - £y?rr +1 

i=l 



2,V, 



y! r? 



Define 



2 

cO 



•yfo-l 



' Sfl o ' 




" V3S' 


0 0 







(70) 
(71) 

(72) 
(73) 



H co = (Yco y 
Hoi = [y c \ y 2 cl 

where SVD of H c0 yields 

H c0 = USV T 

= U«Ej/ 2 Sf V £ (74) 
with R = ran/c(H c0 ). Hence, the reduced-order model is given by 

D 0 = the first N t columns of Y° (75) 

Di = the second N { columns of Y° (76) 

C * U ft Ef (77) 

B„ ~ the first iV,- columns of £^ 1/2 U£ Y 1 (78) 

Bi ~ the second iV,- columns ofE^U^Y 1 (79) 

A ~ E« 1/2 U?H cl V fl S« 1/2 (80) 

In yet another particular embodiment, a method in accordance with the invention may 
be applied to an aerodynamic ROM with only the displacements as the system inputs. This is 
5 possible by applying simultaneously the pulse and double pulse inputs, 



< = S n = 



6? = S" 



1 (n = 0) \ 

0 (n = 1,2,. ..,M) J 

S7 (» = 0) 

"A (» = 1) 

0 (n = 2,3,..., Af) 



(81) 
(82) 
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and get the corresponding responses y^Ydi at the first two time steps: 

Y2 = [y° dl yj 2 ... yj*] (33) 

= [y^ yj, ... y^vj (34) 

For a new SCI, we use 

Vsa = E b oi^ + E b ii^ (35) 

i=l t=l 

d? c/ = Edoitf + £di,r? (86) 

1=1 1 = 1 

where f? are discrete- time derivative of rj\ To be consistent with the use of the 
double pulse defined in (82), r? are obtained by filtering rf via <5f, i.e., 

r? = cont/(r? f #) (87) 



which is equivalent to the backward difference equation, 



rf = U /' (88) 



Subject to the new SCI we sample the system response y" and get 
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